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ABSTRACT 

A finite element formulation is presented for the analysis of large deflection response of composite panels 
subjected to aerodynamic pressure at supersonic flow and high acoustic excitation. The first-order shear 
deformation theory is considered for laminated composite plates, and the von Karman nonlinear strain- 
displacement relations are employed for the analysis of large deflection panel response. The first-order piston 
theory aerodynamics and the simulated Gaussian white noise are employed for the aerodynamic and acoustic loads, 
respectively. The nonlinear equations of motion for an arbitrarily laminated composite panel subjected to a 
combined aerodynamic and acoustic pressures are formulated first in structure node degrees-of-freedom. The 
system equations are then transformed and reduced to a set of coupled nonlinear equations in modal coordinates. 
Modal participation is defined and the in-vacuo modes to be retained in the analysis are based on the modal 
participation values. Numerical results include root mean square values of maximum deflections, deflection and 
strain response time histories, probability distributions, and power spectrum densities. Results showed that 
combined acoustic and aerodynamic loads have to be considered for panel analysis and design at high dynamic 
pressure values. 
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1. INTRODUCTION 

Aircraft and spacecraft skin-panels are subjected simultaneously to high levels acoustic (sonic fatigue) and 
aerodynamic (panel flutter) pressures. 1 ^ Sonic fatigue and panel flutter have been the major design considerations 
for aircraft, spacecraft and missiles since the late nineteen sixties. An excellent review of sonic fatigue technology 
up to 1989 was given by Clarkson. 3 Various types of pressure loads, developments of theoretical methods, and 
comparisons of experimental results with theories and nomographs were given. Recently, Wolfe et al. 4 gave 
reviews in-depth of sonic fatigue design guides, classical and finite element approaches, and identification 
technology. Experimental investigation of nonlinear beams and plates, and multimodal fatigue model were also 
reported. Sonic fatigue design guides have been developed by Rudder and Plumblee 5 for isotropic metallic and by 
Holehouse for graphite-epoxy composite aircraft structures. The design guides, however, were based on the semi- 
empirical test data or the simplified single-mode approach. Vaicaitis et al. have developed a Galerkin-like 
procedure (PDE and modal method) and a time domain Monte Carlo approach for the nonlinear response of 
isotropic 1,2 and composite 7 8 panels to acoustic and thermal loads. 

An excellent survey of nonlinear panel flutter through 1970 was given by Dowell. 9 The vast amount of 
theoretical literature on panel flutter were grouped into four categories based on the linear or nonlinear structure 
theories, and the two aerodynamic theories (quasi-steady first-order piston or full linearized inviscid potential 
flow). The partial differential equations (PDE), the Galerkin’s method and the numerical simulation results showed 
that a minimum of six modes are needed for a converged limit cycle amplitude response. Bismarck-Nasr 10 gave a 
review of the linear panel flutter using the finite element methods. Recently, a review of various analytical methods 
and experimental results of nonlinear panel flutter at supersonic and hypersonic speeds was given by Mei et al. 11 An 
approach for the design of flutter-free surface panels using the quasi-static Ackeret aerodynamic theory was 
documented by Laurenson and McPherson. A exhaustive search of the literature reveals that there is no study of 
nonlinear panel response to combined acoustic and aerodynamic loads. 

When a flight vehicle travels at supersonic speeds, panel flutter caused by aerodynamic pressure is not the only 
form of dynamic instability. The surface panels also experience high frequency random pressure fluctuations (sonic 
fatigue). This paper presents a finite element formulation for the analysis of nonlinear large deflection response of 
composite panels subjected to high acoustic excitation and aerodynamic pressure at supersonic speeds. The first- 
order shear deformation theory is considered for the laminated composite plates. The von Karman nonlinear strain- 
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displacement relations are employed for the large deflection response of the panel. Simulated Gaussian white noise 
and the first order piston theory aerodynamics are employed for the acoustic and aerodynamic loads. The nonlinear 
equations of motion for an arbitrarily laminated composite panel subjected to a combined high acoustic and 
aerodynamic loads are formulated first in the structure node degrees-of-freedom (DOF). The system equations are 
then transformed and reduced to a set of coupled nonlinear equations in modal coordinates. Numerical integration is 
employed to obtain the panel response. Examples are given for an isotropic and a composite panel at various 
combinations of sound pressure level and dynamic pressure. 


2. FORMULATION 

2.1 Equations of Motion in Structure Node DOF 

The inplane strain, curvature and shear strain vectors based on the von Karman large deflection and the first 
order shear deformation theories are given by 
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where u, v and w are the inplane and transverse displacements, respectively. ¥* and 'Ey are the rotations of the 
normal to the midsurface about the y- and x-axes, respectively. The subscripts m and b denote membrane (inplane) 
and bending components, respectively. The constitutive equations for a laminated composite plate are 

U ° ] {cMaM ( 2) 

where [A], [B], [D] and [A s ] are the stretching, bending-stretching coupling, bending, and shear stiffnesses, 
respectively. The quasi-steady first-order piston aerodynamic theory is employed for the aerodynamic pressure at 
high supersonic Mach number (M«> 1 .6). The aerodynamic pressure is given by 9 
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where q a =p a V« 2 /2 is the free stream dynamic pressure, p a the air density, V_ the velocity and /3 = -JmI -1 

Using the Hamilton’s principle and the finite element expressions, the system equations of motion for a 
composite plate subjected to aerodynamic pressure and high acoustic excitation can be expressed as 1314 

^["l{>n}+^-[c]{w,} + ([Kj + [ir»Mi;}={P t (,)} ( 4) 

where a>o - ^£> 110 / pha 4 j 1 1 2 is a reference frequency and p, h and a are the panel density, thickness and length, 

respectively. The value D 110 is the first entry in laminate bending rigidity [D] calculated when all of the fibers of the 
composite layers are aligned in the airflow x-direction. The [M] b , [G] and (P b ) are the system mass matrix, 
aerodynamic damping matrix and load vector due to random acoustic pressure, respectively. The linear and 
nonlinear system stiffness matrices are given by 

fo]= ]+ M* + [Ks ]- [k b \k]^ [KbY 

l K NL]—[KB\KXn ^^mb + \-^B\ + [ K ^Nm] + [^Nb] + [^\~[^]bm\XXn{^BY + [^ m b j (5) 

where the nondimensional dynamic pressure and nondimensional aerodynamic damping are given by 
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where C a = ji(mI - 2\ I is the aerodynamic coefficient and p = p o alph is the mass ratio. For high 
supersonic speeds M 00 »l, Dowell 9 approximated C fl = pf M M . [A a ] is the system aerodynamic influence matrix, 

and [Kl] and [K2] are the first-order and second-order nonlinear stiffness matrices which depend linearly and 
quadratically upon the unknown system displacement vector {W}={{W b },{W m }} T The subscripts B, N m and N b 
denote that the corresponding stiffness matrix is due to the laminate extension-bending stiffness [B], membrane 
force components {N m }=[A]{e m 0 } and {N b }==[B]{K}, respectively, and the subscripts a and s denote aerodynamic 
and shear deformation, respectively. The derivation of Eq. (4) can be referred to references 13 and 14. In the 
absence of acoustic loading, {P b (t)}=0, Eq. (4) reduces to nonlinear panel flutter at supersonic speeds; On the other 
hand, by setting k=g a =0, Eq. (4) describes nonlinear random response of a composite panel subjected to high 
acoustic excitations. According to the authors’ knowledge, this is the first attempt in investigating of nonlinear 
response of panels subjected to a combined acoustic and aerodynamic pressures. 

For a given set of X and C a (or Eq. (4) can be solved by numerical integration in the structure node 

DOF for combined load case. This approach has been carried out for sonic fatigue analysis with simulated random 
loads. ’ It turned out to be computationally costly due to: (i) at each time step, the element nonlinear stiffness 
matrices are evaluated and the system nonlinear stiffness matrix [Knl] is assembled and updated, (ii) the number of 
structure node DOF of { W b } is usually very large, and (iii) the time step of integration should be extremely small. 

An efficient solution procedure is to transform Eq. (4) into the modal coordinates with a modal reduction. This 
approach is presented as follows. 

2.2 Equations of Motion in Modal Coordinates 

Express the panel deflection as a linear combination of some known base functions as 

tybh StfrWWMofo} (7) 

r=l 

where the number of retained linear in-vacuo modes, n, is much smaller than the number of structure node DOF in 
bending, {W b }. The normal mode { <J) r } , which is normalized with the maximum component to unity, and the linear 
natural frequency C0r are obtained from the linear vibration of the system 

k l*r }- (ft + IK, ]- [Kb I* fc 1 [*„ f ) (8) 

A small number of most contributing modes to be retained in the analysis can be determined from the modal 
participation value which is defined as 

n 

Participation of the r mode = RMS(q r )/ RMS(q s ) (9) 

s = 1 

Since matrices [Kl]^, [KIb], [Kl^] and [K2] are all functions of the unknown bending DOF {W b }, they can 
now be expressed as the sum of products of modal coordinates and nonlinear modal stiffness matrices as 

I ] (10) 

r=l [ s = 1 ; 

where the super-indices of those nonlinear modal stiffness matrices denote that they are assembled from the 
corresponding element nonlinear stiffness matrices. Those element nonlinear stiffness matrices are evaluated with 
the corresponding element components {w b } (r) obtained from the known system linear mode {<f) r }. Therefore, the 
nonlinear modal stiffness matrices are constant matrices. The matrix [KInoJ, however, is a linear function of the 
inplane DOF { W m } which consists of two terms as 
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where the two inplane modes corresponding to the r Ih bending mode {<|) r } are given by 

= {♦.L-W:inEk} m 

The nonlinear stiffness matrix [Kl Nm ] can be expressed as the sum of two nonlinear modal stiffness matrices as 

[«».]— Z IrfMjfcF - 1 iq rqs [K2 Nm f rS) (13) 
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The nonlinear modal stiffness matrices [Kl Nn J (r) and [K2 Nn J (rs) are constant matrices and they are assembled and 

evaluated with the known inplane modes {<}) r } m and {<t> rs } m , respectively. Equation (4) is thus transformed to the 
reduced modal coordinates as 

^-Cirl{','} + + it, a. - (Prl+ [*,]+ [*„]){,}= {?«)} <M> 

where the modal matrices are given by 
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and the quadratic and cubic terms are 

[*« f® f Hm. ([«]'"> - [fin, f rs> - [m]« [*£ [*ri]W W,} 

and the modal force is 
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A structural modal damping 2 !&©, ~^[ 1 ) has been added to Eq. (14), and £ is the modal damping ratio which can 

be determined experimentally or from the data base of structures of similar construction. The nonlinear response for 
a given panel at certain dynamic pressure X and damping parameters C, and ^ can be determined from Eq (14) by 
any numerical integration scheme. The advantages in using Eq. (14) are: (i) there is no need to assemble and update 

the quadratic and cubic nonlinear terms since all the nonlinear modal matrices are constant matrices and (ii) the 
number of modal equations, n, is small. 

2.3 Random Surface Pressure 

The input acoustic excitation is assumed to be band-limited Gaussian random noise and uniformly distributed 
over the structural surface. The power spectrum density (PSD) has the form 

s (f) = Po 10 SPL,1 ° 0 <f<f c 
= 0 otherwise 

where p 0 is the reference pressure, p o =2.9xl0‘ 9 psi (0.00002 N/m 2 ), SPL is the sound pressure level in 
decibels and f c is the selected band width. The formulation presented in Eq. (14), however, is not limited 
to stationary Gaussian excitation. It can also handle nonstationary, non-Gaussian random loading which 
the high speed flight vehicles would probably experience. With recorded flight high frequency pressure 
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fluctuations, random panel response can be determined much realistically by numerical integration of Eq. 
(14). 

3. EXAMPLES AND DISCUSSION 

The nonlinear system equations presented in Eq. (4) are general in the sense that they are applicable for 
rectangular 1 15 or triangular 14 finite elements. The finite element employed in the examples is the three-node 
triangular Mindlin (MIN3) plate element with improved transverse shear. 1 The shear correction factor is defined 
as a s =1/(1 +0.5 2^/ 2X„) . The MIN3 element has a total of 15 DOF, 5 at each apex node. The bending node DOF 

r=4.9 i=4 , 9 

{w b } comprise of transverse displacements and normal rotations (w, *F X and 'Ey) and the inplane node DOF {w m } 
comprise of inplane displacements (u and v). Nonlinear response are obtained for a square isotropic plate and a 
rectangular composite plate. An aerodynamic coefficient C a =0.01 and a modal damping ratio £,=0.01, r=l to n are 
used in the examples. 

3.1 Square Isotropic Plate 

A simply supported square aluminum panel under the combined acoustic and aerodynamic pressures is studied 
in detail. The plate is of 12x12x0.040 in. (30.5x30.5x0. lcm) with immovable inplane edge conditions u(0,y) = 
u(a,y)= v(x,0)= v(x,a)= 0. The material properties are E=10Msi (68.89 GPa) and v=0.3. The plate is modeled with a 
12x12 mesh or 288 MIN3 elements. The number of structural node DOF {W b } is 407 for the system equations 
given in Eq. (4). It is well know that 6 modes in the airflow direction are needed for a converged limit cycle 
response for panel flutter and the lowest few symmetrical modes are needed for the uniform input random pressure 
distribution for sonic fatigue. Therefore the modes considered for the combined aerodynamic and acoustic loads are 
(1,1) to (6,1) for the panel flutter and (1,1), (1,3), (3,1) and (3,3) for the sonic fatigue analysis. No modal 
participation calculations are needed for this well studied problem in panel flutter and sonic fatigue. The root mean 
square (RMS) maximum deflections to plate thickness versus the nondimensional dynamic pressure X at SPL of 0, 
100, 1 10 and 120 dB are shown in Fig. 1. 



Fig. 1 RMS maximum deflection of a simply Fig. 2 Comparison of limit cycle amplitude for a 

supported square isotropic plate simply supported square isotropic plate 


The curve with the null acoustic pressure (0 dB SPL) is the conventional panel flutter limit cycle oscillations, 
and the RMS deflections at null dynamic pressure (X=0) are the conventional nonlinear panel response to acoustic 
excitations. The accuracy of the present formulation for panel flutter and sonic fatigue is verified and shown in Fig. 
2 and Table 1 for panel flutter and sonic fatigue, respectively. The Fokker-Planck-Kolmogorov (FPK) Equation 
[18] method is an exact solution to the single degree-of-freedom forced Duffing equation. The finite 
element/equivalent linearization (FE/EL) approach assumes that the equivalent linearizd system obtained after the 
application of equivalent linearization technique is also stationary Gaussian, while the present time domain 
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numerical integration method does not assume that the displacement response is Gaussian, therefore, the present 
method should be more accurate and realistic. 

Table 1. Comparison of RMS of (W„„/h) for a simply supported rectangular 
(15xl2x0.040Jm)jsotroj)icj>latejisingjUfferentinethod^andnumber of modes 

SPL FPK [18] FE/EL [19] Present 
(dB) 1 mode 1 mode 4 modes 4 modes 
90 0.249 0.238 | 0.238 0.266 

100 0.592 0.532 0.533 0.489 

110 1.187 1.030 1.031 1.092 

120 1 2.200 I 1,902 | 1,905 | 2.113 

The maximum deflection is located at the three-quarter length from the leading edge. (3a/4, a/2) for the panel 
flutter, however, it is at the plate center (a/2, a/2) for sonic fatigue. For a fixed SPL, the location of the maximum 
deflection thus moves from the plate center towards the three-quarter length as the dynamic pressure increases. On 
the other hand, for a fixed dynamic pressure, the location of the maximum deflection moves towards the center 
from the three-quarter length by increasing the SPL. The maximum deflection for the combined acoustic and 
aerodynamic loading case, therefore, is not at one fixed location, it can be anywhere between the plate center and 
three-quarter from the leading edge. 

The deflection results shown in Fig. 1 indicate that: (i) the superposition method does not applicable to 
nonlinear system, (ii) for a given SPL, the RMS deflection at X=0 is higher than the deflection at 0<X<X cr , this is 
due to the increase of panel frequency with increasing of X, thus increase the stiffness of the panel, and (iii) the 
RMS deflection at large dynamic pressure (X»X cr ) are always higher than those at X=0. Thus important 
conclusions can be drawn for design and analysis of surface panels at supersonic speeds: (i) for X«X cr , only 

acoustic loading or sonic fatigue has to be considered, and (ii) for X>X„, both acoustic and aerodynamic loads have 
to be considered. 

Representative panel behavior at five loading combinations are presented in Figs. 3-7. Those five loading 
combinations correspond to A to E shown in Fig. 1. The maximum deflection and maximum strain response time 
histones, probability distnbution and PSD for each loading case are presented. Figures 3 and 4 show the random 
response at SPL=100 and 120 dB and X=0 (sonic fatigue; points A and B in Fig. 1), respectively. At the low 100 dB 
SPL, the panel basically experiences a small deflection linear random vibration dominated by the fundamental (1,1) 
mode. The panel motion at the high 120 dB SPL, however, is clearly a large deflection nonlinear random vibration. 
This is demonstrated by the peaks in PSD plots that are broadening and shifting to the higher frequency and by the 
presence of nonzero-mean inplane strain shown in strain plots. 

At X=800 and 0 SPL (panel flutter; point C in Fig. 1), the panel. response shown in Fig. 5 is a large amplitude 
limit cycle motion. The displacement probability density and the PSD of(W rax /h) both describe a periodic motion. 

The maximum strain time history shows clearly the effect of large inplane strain component due to large amplitude 
periodic motions. 

The panel responses at the combined loads of X=800 and SPL=100 and 120 dB (points D and E in Fig. 1) are 
shown in Figs. 6 and 7, respectively. The maximum deflection and strain time histories show the nonlinear large 
deflection vibrations dominated by the fundamental mode and the presence of inplane strain components. 

3.2 Rectangular Composite Plate 

Nonlinear response of composite panels under aerodynamic and acoustic pressures can be determined using 
the present formulation and solution procedure. As an example, a clamped rectangular graphite-epoxy plate of eight 
layers [0/45/-45/90]s is analyzed. The dimensions and material properties of the panel are: 
a =15 in. (38.1 cm) E,=22.5 Msi (155 GPa) 

b = 1 2 in. (30.5 cm) E 2 = 1 . 1 7 Msi (8.07 GPa) 

h=0.048 in. (0. 122 cm) G, 2 =0.66 Msi (4.55 GPa) 

p=0. 1458x1 0' 3 lb-s 2 /in 4 (1550 Kg/m 3 ) G 23 =0.44 Msi (3.03 Gpa) 

v 12 =0.22 
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The inplane edges are immovable and the plate is modeled with a 12x12 mesh. The number of system equations in 
structure node DOF { W b } is of 363. The system equations are reduced to the modal coordinates using the lowest n 
modes in increasing frequency order. The RMS W^/h at 120 dB SPL and A^=800 using different number of modes 
are shown in Table 2. The results show that a 20- or 25-mode model would yield a converged RMS maximum 
deflection. 

To demonstrate the advantage of using modal participation defined in Eq. (9), the participation values for the 
25-mode model are shown in Table 3. By retaining those 13 modes with participation value > 1% in the analysis, 
the RMS Wma X /h is 0.8124 at 120 dB SPL and yU=800. Using the 13 most contributing modes, the RMS (Wmax/h) 
versus nondimensional dynamic pressure X at SPL of 0, 100, 110 and 120 dB are shown in Fig. 8. Similar 
conclusion from the isotropic panel can be drawn for the composite panel, that is, at low dynamic pressure (X,«>l cr ) 
only acoustic loading needs to be considered, and at high dynamic pressure (^>Xc r ) both aerodynamic and acoustic 
loads have to be considered for the design and analysis of surface panels at supersonic flow. Response time history, 
probability distribution and PSD are not repeated for the composite panel. 

Table 2. RMS (W^x/h) for a clamped rectangular graphite-epoxy panel 
at 120 dB SPL and X=800 using different number of modes 


Number of modes, n 

RMS (W m /h) 

i 

0.5557 

2 

0.5845 

6 

0.7814 

9 

0.7798 

16 

0.8279 

20 

0.8110 

25 

0.8183 

Selected 1 3 modes 

0.8124 


Table 3. Modal participation values for a clamped rectangular graphite-epoxy 
panel at 120 dB SPL and X=800 using the lowest 25 modes 


Mode number 

1 

2 

3 4 

5 

6 

7 

8 

9 

10 

11 

12 

Participation, % 

36.72 5.24 19.30 4.25 

4.01 

7.67 

1.71 

0.76 

0.33 

1.54 

4.77 

0.35 

13 

14 

15 

16 17 

18 

19 

20 

21 

22 

23 

24 

25 

0.28 

4.19 

0.38 

0.99 1.38 

2.55 

0.21 

0.27 

0.53 

0.33 

0.54 

0.14 

1.54 


CM 



Fig. 8 RMS maximum deflection of a clamped rectangular graphite-epoxy plate 
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4. CONCLUDING REMARKS 

A finite element time domain modal formulation is presented for the analysis of nonlinear response of 
composite panels subjected to combined acoustic and aerodynamic pressures. The advantage of using modal 
participation for retaining the most contributing modes was demonstrated. For panels at supersonic flow, only 
acoustic excitations (sonic fatigue) are to be considered for \«X CT , and both acoustic and aerodynamic pressures 
have to be considered for X>Xcr. Future extension of the present work includes the combined acoustic, aerodynamic 
and thermal loads. Arbitrary flow directions and curved panels will also be considered. 
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Fig. 3 Random Response of a simply supported square isotropic plate at 100 dB SPL and X=0 
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Fig. 7 Random Response of a simply supported square isotropic plate at 120 dB SPL and X=800 


199 

















